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We investigate the spectrum of the staggered Dirac operator in 4d quenched U(l) lattice gauge theory and 
its relationship to random matrix theory. In the confined as well as in the Coulomb phase the nearest-neighbor 
spacing distribution of the unfolded eigenvalues is well described by the chiral unitary ensemble. The same is true 
for the distribution of the smallest eigenvalue and the microscopic spectral density in the confined phase. The 
physical origin of the chiral condensate in this phase deserves further study. 
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By now it is a well-known fact that the spec- 
trum of the QCD Dirac operator 
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in a chiral basis 



is related to universality classes of random ma- 
trix theory (RMT), i.e. determined by the global 
symmetries of the QCD partition function [Q. 
In RMT the matrix T in Eq. (|l|) is replaced by 
a random matrix with appropriate symmetries, 
generating the chiral orthogonal (chOE) , unitary 
(chUE), and symplectic (chSE) ensemble, respec- 
tively 1). For SU(2) and SU(3) gauge groups 
numerous results exist confirming the expected 
relations 

We have investigated 4d U(l) gauge theory de- 
scribed by the action 5{J7;} = ^p{^—cos9p) with 
Ui = Ux,f, = exp(z6'j;,^) and 9p = 9x,f_c + Ox+p.^i^ - 
Ox+o,^ - Ox^u (J^ m) • At I3c « 1.01 U(l) gauge 
theory undergoes a phase transition between a 
confinement phase with mass gap and monopole 
excitations for (3 < (3c and the Coulomb phase 
which exhibits a massless photon Q for (3 > (3c- 
The question of the order of this phase transi- 
tion, and hence the issue of a continuum limit 
/3 — > /3c — of the massive phase, has remained a 
subject of endless debate, see for instance H] and 
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references therein. For (3 > (3c & critical line of 
continuum theories ought to exist as the photon 
is massless for all these /3- values. 

We are interested in the relationship between 
U(l) gauge theory and RMT across this phase 
transition The Bohigas-Giannoni-Schmit j^] 
conjecture states that quantum systems whose 
classical counterparts are chaotic have spectral 
fluctuation properties (measured, e.g., by the 
nearest-neighbor spacing distribution P{s) of un- 
folded eigenvalues) given by RMT, whereas sys- 
tems whose classical counterparts are integrable 
obey a Poisson distribution, P(s) — exp(— s). As 
the Dirac operator is a quantum-mechanical ob- 
ject without classical counterpart, the meaning of 
the Bohigas-Giannoni-Schmit conjecture for lat- 
tice gauge theory is somewhat unclear. Neverthe- 
less, as for SU(2) and SU(3) gauge groups, we ex- 
pect the confined phase to be described by RMT, 
whereas free fermions are known to yield the Pois- 
son distribution. The question arose whether the 
Coulomb phase will be described by RMT or by 
the Poisson distribution, with apparently no clear 
theoretical arguments in favor of either scenario. 

In Ref. Q some of the authors have resolved 
this question by numerical analysis. We gener- 
ated twenty (or more) gauge configurations at 
each of the following parameter values: 8"^ x 4 
lattice at /3 = 0, 0.9, 0.95, 1, 1.05, 1.1, 1.5 and 
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Figure 1. Nearest neighbor spacing distribution P{s) on an S'^ x 6 lattice in the confined phase (left plot) 
and in the Coulomb phase (central plot), and for the free Dirac operator on a 53 x 47 x 43 x 41 lattice 
(right plot). The theoretical curves are the chUE result, P{s) = 32(s/7r)^ exp(— 4s^/7r), and the Poisson 
distribution, P{s) = exp(— s). 



8^ X 6 lattice at /3 = 0.9, 1.1, 1.5. The nearest- 
neighbor spacing distributions for the 8^ x 6 lat- 
tice at /3 = 0.9 (confined phase) and at /3 = 1.1 
(Coulomb phase), averaged over 20 independent 
configuration, are depicted in Fig. |l|. Both are 
well described by the chUE. In contrast, the right 
plot in Fig. |l| shows that free fermions are de- 
scribed by the Poisson distribution. The large 
prime numbers for the lattice size are needed to 
remove the degeneracies of the spectrum. 

We have continued the above investigation with 
a study of the distribution of small eigenvalues 
in the confined phase. The Banks-Casher for- 
mula [|j relates the eigenvalue density p(A) at 
A = to the chiral condensate, 

S = = lim lim TTp{0)/V . (2) 
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The microscopic spectral density 

should be given by the result for the chUE 
of RMT 0]. This function also generates the 
Leutwyler-Smilga sum rules j|] . 

To study the smallest eigenvalues, spectral av- 
eraging is not possible, and one has to produce 
large numbers of configurations. Our present re- 
sults are for (3 = 0.9 in the confined phase with 
10000 configurations on a 4**, 10000 configuration 
on a 6^, and 1106 configurations on an 8"^ x 6 lat- 
tice. The upper plot in Fig. ^ exhibits the distri- 
bution P(Amin) of the smallest eigenvalue Amin in 



comparison with the prediction of the (quenched) 
chUE of RMT for topological charge v — 0, 



-P(Amin) 



exp 



(FEA„,in)^ 



(4) 



The agreement is excellent for all lattices. For 
the chiral condensate we obtain S « 0.35 by ex- 
trapolating the histogram for p(A) to A = and 
using Eq. (^). (With staggered fermions on a fi- 
nite lattice one always has p(0) = 0, but within 
reasonable limits the extrapolated value is inde- 
pendent of the choice of the bin size.) Since the 
average value of Amin goes like V~^, (Amin) de- 
creases with increasing lattice size. In the lower 
plot of Fig. I the same comparison with RMT is 
done for the microscopic spectral density (^) up 
to z = 10, and the agreement is again quite satis- 
factory. Here, the analytical RMT result for the 
(quenched) chUE and ~ is given by 



(5) 



The quasi-zero modes which are responsible for 
the chiral condensate S w 0.35 build up when we 
cross from the Coulomb into the confined phase. 
For our 8^ x 6 lattice, Fig. | compares on identical 
scales densities of the small eigenvalues at /? = 0.9 
(left plot) and at (3 = 1.1 (right plot), averaged 
over 20 configurations. The quasi-zero modes in 
the left plot are responsible for the non-zero chi- 
ral condensate S > via Eq. (||), whereas no 
such quasi-zero modes are found in the Coulomb 
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Figure 2. Distribution P(Ainin) (upper plot) and 
microscopic spectral density Psiz) (lower plot) 
from our 6"* lattice data in comparison with the 
predictions of the chUE of RMT (dashed lines), 
see Eqs. (|) and (|). 

phase. This is as expected. However, it may be 
worthwhile to understand the physical origin of 
the U(l) quasi-zero modes in more detail. For 
4d SU(2) and SU(3) gauge theories a general in- 
terpretation is to link them, and hence the chiral 
condensate, to the existence of instantons. As 
there are no instantons in 4d U(l) gauge theory, 
one needs another explanation, and it may be in- 
teresting to study similarities and differences to 
the 4d SU(2) and SU(3) situations. An analogous 
case exists in 3d QCD po|] . 

In conclusion, the nearest-neighbor spacing dis- 
tribution of 4d U(l) quenched lattice gauge the- 
ory is described by the chUE in both the confine- 
ment and the Coulomb phase. In the confinement 
phase we also find that the P(Ainin) distribution 
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Figure 3. Density p{X) of small eigenvalues for 
the 8^ X 6 lattice at /3 = 0.9 (left plot) and at 
/3 — l.l (right plot). A non-zero chiral condensate 
is supported in the confinement phase. 



and the microscopic spectral density (^ are de- 
scribed by the chUE. A better physical under- 
standing of the origin of the quasi-zero modes, 
which are responsible for the non-zero chiral con- 
densate, is desirable. 
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